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Abstract 

This  paper  describes  results  on  open-  and  closed-loop  controlled  growth  of  epitaxial  GaP  / 
Ga]_xInxP  heterostructures  on  Si(001)  substrates.  The  layers  are  grown  in  a  low  pressure 
pulsed  chemical  beam  epitaxy  (PCBE)  reactor  utilizing  real  time  optical  p-polarized 
reflectance  (PRS)  probing.  The  results  of  the  implemented  closed  loop  controlled  growth 
favorably  compare  to  the  films  grown  using  pre-designed  source  injection  profiles  based  on 
an  experimental  data  base. 
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INTRODUCTION 


The  development  of  surface- sensitive  optical  characterization  techniques  of  thin  film 
growth  opened  up  the  possibilities  to  obtain  detailed  information  in  real  time  that  is  highly 
relevant  to  the  growth  kinetics.  The  significance  of  these  techniques  for  growth  monitoring  is 
that  they  move  the  observation  point  close  to  the  surface  of  the  film  where  the  growth  occurs. 
These  developments  make  it  feasible  to  improve  thin  film  deposition  applying  closed  loop 
control  techniques.  Areas  of  interest  include  fabrication  of  advanced  nanostructure  devices, 
improved  densities  of  integrated  electronic  devices. 

In  this  work  we  describe  our  monitoring  and  control  method  for  thickness/composition 
control  of  S i (00 1  )/Ga P/Ga j  _x  I  n  x  P  growth  in  low  pressure  PCBE  system.  The  films  grown 
with  our  control  algorithm  have  been  analyzed  by  ex-situ  methods  such  as  x-ray  diffraction 
(XRD)  and  Secondary  Ion  Mass  Spectroscopy  (SIMS).  The  test  results  show  the  superiority 
of  the  closed  loop  controlled  growth  over  open  loop  control  designed  using  our  data  base, 
which  was  built  up  from  accumulated  growth  data.  In  our  experiment  we  use  p-Polarized 
Reflectance  Spectroscopy  (PRS)  as  the  primary  probing  method  of  our  thin  film  growth 
process.  PRS  is  based  on  the  changes  in  reflectivity  R  of  the  p-polarized  laser  light  shone  on 
a  dynamically  changing  stack  of  dielectric  media  and  its  sensitivity  for  the  growth  kinetics  is 
based  on  choosing  the  incident  angle  to  be  the  pseudo-Brewster  angle  (to  Si  and  GaP 
substrates  in  our  case  ). 

In  earlier  papers1-3  we  have  reported  on  our  approach  of  modeling  the  surface  kinetics  of 
epitaxial  GaP  and  Gaj_xInxP  growth  on  Si.  For  both  the  GaP  and  Gaj_xInxP  growth  stages 
we  introduced  reduced  order  surface  kinetics  (ROSK)  models  to  represent  the  essential 
chemical  processes  in  the  surface  reaction  layer  (SRL).  In  this  paper  we  utilize  our  model  to 
monitor  and  control  the  growth  rate  and  composition  of  the  growing  film  as  follows.  First, 
for  the  modeling  of  the  PRS  reflectance  measurement  we  use  Fresnel's  equation  and  a  virtual 
interface  method,  introduced  by  D.E  Aspnes4’5,  for  the  multi-layer  stack  of  GalnP  where  the 
change  of  composition  and  thickness  of  the  growing  layer  is  determined  by  our  model 
dynamics  and  the  flow  rates  are  entered  as  input  variables.  We  formulate  the  control  of  thin 
film  growth  as  an  optimal  control  problem.  Second,  we  use  a  nonlinear  filtering  algorithm6 
to  estimate  growth  rate  and  composition  of  growing  Gaj_xInxP  film  and  based  on  these 
estimates  then  determine  optimal  flow  rates  of  our  source  vapors  to  achieve  the  desired 
composition  and  growth  per  cycle  in  real  time. 

Due  to  the  limited  space  for  a  detailed  description  of  the  experimental  setup  we  refer  to 
previous  publications7-9. 

Figure  1: 

Schematic 
representation  of  a 
precursor  pulsing 
sequence  used  in  the 
growth  of  the  ternary 
compound 
semiconductor  Ga^ 
xInxP  grown  via  the 
organometallic 
precursors  TBP, 

TEG  and  TMI. 
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DATA  BASE  ANALYSIS/OPEN  LOOP  CONTROL  DESIGN 


In  the  PCBE  growth  of  GalnP  film  layers,  the  quality  of  the  grown  film  strongly  depends 
on  ambient  pressure,  temperature,  the  flow  rate  and  the  pulse  timing  of  the  three  source  gases. 
We  adjusted  appropriately  those  by  selecting  optimal  background  pressure,  temperature  and 
the  injection  rate  of  TBP,  as  well  as  pulsing  profile  over  the  various  growth  conditions.  Thus 
the  control  authority  we  used  are  the  flow  rates  of  TEG  and  TMI  source  gases.  A  series  of 
test  runs  was  conducted  to  obtain  data  base  on  the  effect  of  the  changes  in  the  TM  I/TEG 
injection  ratio  on  the  growth  of  GaxInj.xP  film.  To  identify  the  effects  on  the  growth  rate  and 
composition  we  have  grown  GalnP  films  using  set  flows  for  TEG  and  TMI,  i.e.,  u^mi  /  UTEG 
was  set  to  a  constant  value  that  varied  from  0  to  1.2  (with  fixed  nominal  flow  u^eq  =  0.545 
seem)  on  top  of  an  initial  GaP  film.  We  analyze  the  corresponding  the  PRS  data  and  XRD 
measurements  to  establish  the  functional  relationships  x  =  <b(y),  gr  =  \)/(y) ,  where  y 
denotes  TM  I/TEG  flow  ratio,  x  composition  and  gr  the  Gai_xInxP  film's  growth  rate  as 
depicted  in  Figure  3. 

We  utilize  this  data  base  to  perform  the  open-loop  control  synthesis.  That  is,  given  a 
discretization  zk  with  the  thickness  step  size  Az  and  a  corresponding  average  desired 
composition  sequence  xk  we  select  the  flow  ratio  to  be  yk  =  d>~1(xk)  and  the  corresponding 
duration  tk  =  Az  /  vP"i(yk).  Thus  we  can  design  predetermined  flow  rates  of  TMI  and  TEG  to 
achieve  the  desired  composition/growth  profile. 


Figure  2: 

Growth  monitored 
by  PRS  during 
heteroepitaxial 
Gaj  xInxP/GaP  on 
Si(001). 

Evolution  of  the  PR 
signals  .  The  insets 
show  the  fine  structure 
response  at  two 
different  positions  with 
different  TMI: TEG 
flow  ratios  and  different 
PR  responses  to  it.  The 
control  model 
introduced  below  will 
demonstrate  how  the 
optical  PR  response  is 
linked  to  composition 
and  growth  rate 
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ratio  of  TMI:TEG  control  signal 
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Figure  3: 
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REAL-TIME  MONITORING/CONTROL 


In  this  section  we  describe  our  approach  to  utilize  the  real  time  optical  observations  and 
apply  feedback  control  methodology  for  controlling  GalnP  film  growth.  Consider  the  four 
layer  stack  composed  of  ambient  /  surface-reaction  layer  /  film  /  substrate.  We  assume  an 
effective  media  with  the  homogeneous  dielectric  function  £|  in  the  surface  reaction  layer. 
Let  us  denote  the  four  media  by  the  indices  labeled  from  the  ambient  to  substrate.  The 
reflection  coefficient  rn_\  n  from  the  (n-1)  -st  layer  to  n-th  media  is  given  by 


en  i 

/  £n-i  -  eo  sin2(p  -  £n-i  V 

.  2 

8n  -  £q  sin  cp 

£n  V 

8n-i  -  £0  sin2cp  +  en_!  ^ 

1  8n  -  £0  sin2cp 

(1) 


where  £n  is  the  dielectric  constant  of  the  n-th  media.  The  phase  factor  <f>n  for  the  n-th 
media  is  given  by 


2  7i  d 


£q  sin  cp 


(2) 


where  dn  is  the  thickness  of  the  n-th  media. 

For  the  case  of  multi-later  stack  of  films  by  applying  the  theory  of  the  virtual  interface 
method  by  Aspnes10’n,  the  reflectance  amplitude  r  of  the  p-polarized  light  is  then  given  by 


r01  ~ 


-2  i  <t>i 


i  *  -2  i 

1  +  r01  r  e  1 


with  r : 


r12  -  rke 


-2  i  <t>- 
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1  +  r12  rke  2 


(3) 


The  virtual  reflection  index  is  updated  by 


rk  = 


>k,k-l  -  >k-l  e 


-2  i  <I>n 


1  +  %k-l  T-1  e 


-2  i  <t>2 


with  rk  =  \  e  1  0k  , 


(4) 


at  the  end  of  cycle,  where  defines  the  phase  factor.  Based  on  the  phase  factor  we  estimate 
the  thickness  of  the  grown  layers.  For  each  homogeneous  layer  we  have  the  estimate  of  the 
thickness  d2  by 
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do  — 


x 


4  7T-J 


e2  -£osin  <?() 


(  ^cud  '  ^begin)  > 


(5) 


where  the  9enci ,  9begin  is  the  phase  factor  at  the  end  and  beginning  of  the  layer.  Similarly  the 
growth  grk  per  each  cycle  k  is  given  by 


srk  =  - 


4  7t  -J 


e2  -  e0  sin  <p0 


(  &k  '  ^k-l)  • 


(6) 


We  use  the  nonlinear  filtering  algorithm6  for  estimating  the  state  consisting  of  the  virtual 
reflection  index  rk  =  eXl+1X2,  the  film  dielectric  constant  e2  =  x3  +  ix4,  and  growth  per 
cycle  X5  in  real  time.  In  turn  the  thickness  of  the  specific  compound  is  estimated  by  (5).  The 
growth  ratio  of  GaP  and  InP  for  each  cycle  determined  by  (6)  provides  a  composition 
estimate.  Let  y^  denote  the  PRS  signal  at  the  end  of  the  k-th  cycle.  Then  the  filtering  problem 
is  to  estimate  the  signal  process  xk  defined  by 


f  k\ 

7i(*V 

X I 

k 

x2 

fi(xk) 

k 

k- 1 

x3 

— 

k 

k- 1 

x4 

X4 

k 
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\x5) 

l  xs  ; 
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based  on  the  observation  process  yk  =  h(xK )  +vk.  Here  we  assumed  that 


\rk,  k-i 


is 

be 


sufficiently  small  and  used  rk  =  rk_{  c"  -  for  updating  the  virtual  index  r^.  If  we  let  Jr 
the  growth  ratio  of  GaP  or  InP  to  each  nominal  flow  rate,  then  the  functions/},  /}  and  h  are 
defined  by 


/1  +  i  /2  -  *1  +  i  x2  +  2  i  cp2;  h 


ro2  +  rv 
1  +  r02rv 


(8) 


where  cl2  =  Jr  x5.  We  assume  that  noise  processes  wk,  vk  are  independent  (identically 
distributed)  Gaussian  random  variables  with  mean  zero  and  covariance  Q  and  R ,  respectively. 
The  growth  of  GaP  and  InP  is  determined  in  terms  of  riQaP  and  nInP  which  are  given  by 


d  n 


GaP 


dt 


=  kA 


d  n 


np  nGa  > 


InP 


dt 


=  k5  np  nIn 


(9) 


where  nP,  riQa,  and  nIn  denote  the  concentration  of  surface  active  phosphorous,  gallium  and 
indium,  respectively.  We  consider  the  model  for  the  concentration  change  of  active  Ga  in  the 
SRL  by 

nGa  =  UTEG  $GaP  '  nGaP  ’  (10) 


where  SqciP  is  a  pre-determined  constant.  Integrating  the  first  equation  in  (9)  we  obtain 
nGap(tk+l  )  =  e  ( nGap(tk  )  '  $GaP  UTEG  )  +  ^GaP  « TEG 


(11) 
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r  tk  +1 

where  tk  is  the  starting  time  of  the  k- th  cycle  and  C  =  k4  \  nJt).  The  rate  constant  k 
varies  and  we  estimate  it  in  real  time.  We  use  our  filtering  algorithm  to  estimate  the 
concentration  nk  of  nGap  and  the  accumulated  rate  constant  Ck  for  the  k- th  GaP  cycle  based 
on 


nk-l  +  Ck-1  ( UTEG  '  nk-l) 

kQ  j 

\  Ck-i  ) 

+ 


wk 


(12) 


with  measurement  grk  =  VGaP  nk  +  v  -  .  Here  grk  is  the  growth  rate  of  k- th  GaP  cycle, 
determined  by  Equation  6.  The  growth  of  the  InP  is  modeled  analogously. 

We  determine  the  input  flow  rates  u'}EG  and  uPM1  by  performing 

2 


min  k 

UTEC 


+  Zk)nGaP  -  grd 


2 

o\  k  k—1 

grd 

+  / 

■>  \  UTEG  ~  Up  EG 

+ 

2 

min  k 

UTM\ 

nInP  Zk 

k  k-1 

nGaP  1  _  Zk 

Up  Ml  -  UjMI 

(13) 


Figure  4: 

Control  of 
heteroepitaxial  Ga( 
xInxP  growth:  The 
control  design  consists 
of  three  elements: 

(1)  ROSKM 
described  by  f; 

(2)  Filter  gains  Gj(t) 
based  on 
Nonlinear- 
filtering 
techniques  and 

(3)  Feedback  law  K 
based  on 
Dynamical 
programming. 

subject  to 


c  x 


Virtual  Reactor 


< _ 


nGaP  ~  e  TCG  (nG(iP  '  ^GaP  uTEg)+  ^GaP  UTEG 
nInP  =  e  C™‘  (n InP  ’  S]np  UpMI  S[np  UpM, 


(14) 


respectively.  CPEG  and  CPMI  are  the  current  estimates  of  C  for  GaP  and  InP  cycle,  and  zk  is 
the  desired  composition  at  the  k  cycle.  That  is,  we  control  the  growth  rate  by  uTEG  and  then 
by  uTM1  the  composition  for  each  cycle. 
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FINDINGS:  OPEN  LOOP  AND  CLOSED  LOOP  RESULTS 

We  conducted  a  number  of  tests  to  compare  open  loop  and  closed  loop  control 
performance  in  the  GalnP  growth.  We  present  the  test  results  for  specified  thickness- 
composition  profiles,  depicted  in  Figures  5  and  6,  for  single  -  and  for  multiple  parabolic 
graded  GalxInxP  heterostructures. 

The  specifications  include  constant,  linearly  and  parabolically  graded  composition 
segments  in  terms  of  the  film  thickness.  We  varied  the  thickness  of  the  parabolically  graded 

o 

quantum  wells  between  the  200  -  1000  A  range.  Direct  analysis  of  film  thickness  and 
composition  of  the  films  was  conducted  using  SIMS  measurements.  The  calibration  of  the 
SIMS  data  has  been  made  using  constant  composition  samples  with  compositions  measured 
by  XRD,  as  well  as  using  a  linear  estimate  for  the  sputtering  rate  throughout  the  composition 
range.  In  tracking  the  prescribed  composition  profile  the  feedback  control  clearly  proved  to 
be  superior. 


Figure  5: 

SIMS  analysis  of 
closed  loop  (top)  vs. 
open  loop  (bottom) 
controlled  growth 
results  plotted  against 
growth/composition 
specifications. 
Specifications  include 
a  600A,  wide 
parabolic  composition 
graded  structure. 
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Figure  6: 


SIMS  analysis  of 
closed  loop  (top)  vs. 
open  loop  (bottom) 
controlled  growth 
results.  Specifications 
include  a  repetition  of 
200  A,  wide  parabolic 
composition  graded 
structures. 


Figure  7 : 

Flow  control  signals 
expressed  in 
(TMLTEG)  ratio  for 
open-  and  closed- 
loop  controlled 
parabolically  graded 
GalnP 

heterostructures, 
with  200  A  width. 
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